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Abstract. We study the gluon content of a large nucleus (i) in the semi-classical McLerran-Venugopalan 
model and (ii) in the high energy limit as given by the quantum evolution of the Color Glass Condensate. 
We give a simple and qualitative description of the Cronin effect and high-p T suppression in proton-nucleus 
collisions. 
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1 Introduction 

At high energies and/or for large atomic number A, the 
wavefunction of a hadron is expected to be dominated 
by a high density gluonic system. Gluons having occu- 
pation numbers tp of order l/a s , which is the maximal 
density allowed by their mutual interactions, overlap in 
phase space and saturate [1,2]. A strong classical field 
is asso ciate d with the wavefunction and assumes a value 
A ~ V a^a ~ yJTp r~j 1/g at saturation. At the same time 
scattering amplitudes become of order 1 and unitarity 
limits are reached [3]. The problem can be attacked by 
weak coupling methods, since the non-linear phenomena 
"push" gluons to occupy higher momenta, and the satu- 
ration momentum Q s , which is defined as the scale where 
f{Qs) ~ l/ a sj is a hard scale increasing as a power of 
energy in the small Bj or ken- x limit. 

Presumably one of the most complete and modern ap- 
proaches to saturation is the effective theory of the Color 
Glass Condensate (CGC) [4, 5]. Fast moving partons have 
a large lifetime due to time dilation and act as "frozen" 
sources p for the virtual emission of softer gluons. One 
solves the classical Yang-Mills equations to obtain the 
color field A{p) and then an observable O(A) is deter- 
mined by averaging over the possible color sources, with 
a probability distribution WV[p]. Increasing the rapidity 
Y = ln(l/a;), more gluons need to be included in the 
source, and a resummation of a s Y enhanced terms in 
the presence of a background field leads to a functional 
Renormalization Group Equation (RGE) for Wy [p] [4-7] . 
This RGE gives an infinite hierarchy of non-linear cou- 
pled equations, the so-called Balitsky equations [8]. The 
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first one describes the evolution of the scattering ampli- 
tude (T xy ) Y of a color dipole (a;, y) off the CGC and reads 
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with a s = a s N c /ir. The first three terms correspond to the 
BFKL equation [9] in coordinate space [10], while the last 
one accounts for unitarization effects. Eq. (1) can be closed 
by a mean field approximation, that is by allowing the last 
term to factorize 1 , something which should be reasonable 
assuming that the target is a large (A ^> 1) nucleus [15]. 



2 Classical Saturation 

Classical saturation, where there is no small- a; evolution, 
can be realized only in a large nucleus. The A x N c va- 
lence quarks are the sources for the emission of gluons, 
and in the McLerran-Venugopalan (MV) model [2] they 
are assumed to be uncorrelated for transverse separations 
Ax < ^qcd' so that the probability distribution is given 
by the Gaussian [2, 16] 



Wmv[p] oc exp 



l/A 



2 p a {x)p a (x) 



/4 



(2) 



where p\ = 2a s A/R\ ~ A 1 / 3 A 2 is the color charge den- 
sity squared, with Ra the nuclear radius. Even though 
the sources are uncorrelated, the created field A is ob- 
tained from a non-linear equation. Thus, starting from its 
canonical definition, the gluon occupation number tpA is 

1 However, this factorization is not valid in the region where 
the amplitude is very small, namely when T < a 2 s [11-14]. 
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Fig. 1. The gluon occupation number in the MV model. Thick 
(black), solid (blue) and dashed (green) lines show the total, 
CGC and BS quantities respectively. 
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Fig. 2. The integrated gluon distribution in the MV model. 
Thick (black), solid (blue) and dashed (green) lines show the 
total, CGC and BS quantities respectively. 



not the one that we would obtain from a simple super- 
position of the sources. When A > 1, the density can 
be high, sets the magnitude of the saturation scale as 
Q 2 (A) re A 2 A 1 ' 3 \nA > A 2 , and we find [17] 

V A = — r(0,z) + ^ wist (z), z^k 2 /Q 2 s (A), (3) 

with k the transverse gluon momentum and k its magni- 
tude. Here, T is the incomplete Gamma function, while 
the explicit expression for ip^ 1 ^ can be found in [17]. 

The first term, enhanced by l/a s , dominates for all 
z < 1, as shown in Fig. 1. We interpret this compact dis- 
tribution, falling exponentially at large z, as the occupa- 
tion number in the CGC phase. The twist term 2 contains 
the bremsstrahlung spectrum (BS) ipss ~ 1/ z , and is im- 
portant for the large- z behavior, while it remains finite as 
z — > 0. Due to the lack of correlations among the valence 
quarks, a sum rule exists [2, 17-19] 

/ dz [<p A (z) - VbsM ~ >0; (4) 

the integrated distribution is obtained by "summing" over 
the nucleons when Q 2 = ZQ 2 S (A) > Q 2 S (A) (see Fig. 2). 
Thus, the effect of the repulsive interactions in the nucleus 
is just a redistribution of the gluons in momenta. The 
spectra ifA and <^bs become equal at a scale Q C (A) such 
that A 2 < Q 2 C (A) ps a s Q 2 s (A) < Q 2 S (A) and "infrared" 
gluons in excess in the BS spectrum are located at k ~ 
Q S (A) in the MV model. Therefore, the MV spectrum is 
enhanced around the saturation scale, as shown in Fig. 1. 

As an immediate consequence, let us consider the Cronin 
ratio 

K-pA = -rrjz = = z ip A , (5) 

2 The coefficient in front of this term is {a s \n[Q 2 (A)/A 2 ]}~ 1 , 
which is assumed to be equal to one. In fact, in a running 
coupling treatment of the problem, this identification becomes 
natural [17]. 



with tfp the spectrum of the proton, when this is obtained 
from a simple superposition of the gluons emitted by its 
valence quarks. It behaves as 

• "KpA < 1 if z < 1 

• TZpA ~ 0(1/ a.) > 1 if z~l (6) 

• "KpA -> 1 + if z > 1. 

The ratio, shown in Fig. 3, has a maximum at z m = 
0.435+O(a s ) [17]. The maximal value K^g* = 0.281/a s + 
O(const) corresponds to a pronounced peak [17,18,20] 
originating from the compact nature of the nuclear wave- 
function at saturation. This value increases with A (since 
l/a s = \nQ 2 s (A)/A 2 ). 

3 Quantum Saturation 

Now consider the evolution of a hadron to higher energies. 
Its wavefunction contains more and more soft gluons, due 
to the a s Y increase in the available longitudinal phase 

RpA 

2 / ^ 



1.5 



1 



0.5 



1 2 3 4 5 

Fig. 3. The Cronin ratio in the MV model. Thick (black), solid 
(blue) and dotted (red) lines show the total, CGC and twist 
contributions respectively. 
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space, and correlations among the color sources are in- 
duced. The gluon occupation number may be obtained 
from the (averaged over impact parameter) dipole-hadron 
scattering amplitude as 



<p(k,Y) = — ( ^ exp(-ik-r)T(r,Y), 
a s J irr* 



(7) 



where r is the dipole size. In general, one is not able to 
solve Eq. (1) analytically. Only a "piecewise" expression 
for a s Y > 1 is known, and when translated to <p it reads 



<p(k,Y) = { 




_ J — -?£ ln^ + Zi 



Q s 



if k < Q s 

if k > Q s (8) 

if k > Q s , 



where the dominant behavior of the saturation momen- 
turn is Q 2 s (Y)=#Q 2 s (0)exp[a s x(7s)Y/j s ], with x ( 7 ) the 
eigenvalue of the BFKL equation and j s — 0.628 the asso- 
ciated anomalous dimension [1,21-24]. In Eq. (8), A is an 
undetermined constant and Io is a modified Bessel func- 
tion of the first kind. From the first two pieces in this 
equation, it is obvious that the solution exhibits geomet- 
rical scaling [21-27] below and in a certain wide region 
above Q s ', it depends on k and Y only through the com- 
bination k 2 /Q 2 (Y). It is instructive to do a first step in 
the non-linear evolution, valid so long as Y <C To 
the order of accuracy and for momenta k < Q S (A, Y), it is 
enough to evolve only the compact piece in Eq. (3) which 
will add a correction of order Y. This correction contains 
power-law tails which are generated from the tails of the 
evolution kernel. It is clear that, when Y > X/a s all the 
components will be "mixed" and, unlike the classical case, 
in the quantum case there is no compact distribution for 
k < Qs(Y) and no parametric separation between the so- 
lutions above and below Q S {Y), as can be seen in Eq. (8). 

The analysis of the Cronin ratio is not trivial since 
we do not know the solution in the whole k-Y plane. Fur- 
thermore, given a "point" in this plane, the proton and the 
nucleus can be in different phases, e.g. the nucleus could 
be saturated while the proton is still dilute. However, we 
can understand the generic important features. 
• The proton is "less saturated" than the nucleus, since the 
initial proton scale ~ A 2 is much smaller than the initial 
nuclear one Q 2 (A), and therefore the available transverse 
space for the proton is larger. Thus, the proton evolves 
faster than the nucleus and the ratio 1Z p a decreases. For 
example, along the particular line k = Q S (A, Y), one has 



dY 



< 



n 



pA 



(a s A- 1/3 y-^ 



(9) 



• For fixed Y and for extremely high momenta both sys- 
tems are dilute, described by the solution in the double 
logarithmic approximation (the last piece in Eq. (8)), and 




Fig. 4. From top to bottom, the Cronin ratio for Y = 
0, 1/2, 2, below and near the saturation scale. Solid (black) 
lines correspond to an evolved nuclear wavefunction and dotted 
(red) to an unevolved one. 



the above mentioned difference in transverse space is now 
unimportant. The ratio approaches 1 from below, namely 



dTZ 



pA 



dfc 2 



fc 2 >Q=(A,r) 



> & K pA 



A- 2 - 



1 



(10) 



• The sum rule breaks down for any Y > because of the 
correlations induced among the sources. The peak remains 
for a while, since the sum rule is a sufficient but not a 
necessary condition for the existence of a peak. 

• One can follow analytically the evolution of the peak, 
as shown in Fig. 4, until it becomes of order 1, since this 
happens very fast due (again) to the large separation be- 
tween the scales A 2 and Q 2 (A). The nuclear wavefunction 
is almost unevolved, while the proton is still dilute. One 
finds 



13 max 



0(1) when Y = (l/4)ln 2 (l/a s ) < l/a s . (11) 



• Even though smaller than 1, a peak persists under fur- 
ther evolution until Y ~ l/a s , when the power-law tails 
will have "washed-out" the compact piece in the nuclear 
wavefunction; the peak flattens out due to the nuclear evo- 
lution and the ratio becomes a monotonic function of k 2 , 
that is 

d"K p A 



dk 2 



> when Y > l/a s 



(12) 



These features of saturation and the Cronin ratio [17] ex- 
tend previous discussions [19,28], agree with the results 
obtained in numerical solutions [29], and remain qualita- 
tively unaltered under a running coupling treatment [17]. 



4 Epilogue 

Saturation phenomena can play a significant role in de- 
termining the produced particle spectra in high-energy 
heavy ion collisions. In d-Au collisions at RHIC at BNL, 
final state interactions are not important, and with p T and 
rj the transverse momentum and the (pseudo) rapidity of 



4 



D.N. Triantafyllopoulos: From Classical to Quantum Saturation in the Nuclear Gluon Distribution 



1.5 - 



i i 



0.5 - 



r)=0 



h + +h' 



T!=l 



h + +h' 
2 



r\ = 2.2 



ri= 3.2 



1 2 3 4 5 6 

p T [GeV/c] 



2 3 4 5 6 

p T [GeV/c] 



2 3 4 5 6 

p T [GeV/c] 



2 3 4 5 6 

p T [GeV/c] 



Fig. 5. Tie BRAHMS data [30]: Nuclear modification factor for charged hadrons. 



the produced particle, one probes the nuclear wavefunc- 
tion at a value xau — 2|p T | exp[— ry]/y / s NN , where ^/s NN is 
the center of mass energy per nucleon. At current energies 
one expects to reveal classical saturation properties in the 
mid-rapidity region, while quantum saturation should be 
realized in the forward one. Indeed, the CGC predictions 
(and postdictions) seem to be in reasonable qualitative 
agreement with the RHIC data [30] shown in Fig. 5. One 
should keep in mind that the gluon occupation number 
and the particular ratio we studied are not directly mea- 
surable quantities, and therefore any conclusion drawn at 
the quantitative level might be misleading. Nevertheless, 
more "refined" quantities like, for example, the gluon pro- 
duction [19, 29, 31] (and the corresponding ratio) and even 
the charged hadron production [32], are directly related to 
the gluon occupation number and share the same features 
as those presented in the previous sections. It could be 
very well the case that the data for the nuclear modifica- 
tion factor shown in Fig. 5 correspond to a manifestation 
of saturation in the nuclear wavcfunction. 
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